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ABSTRACT 


Magsat data has been analyzed as a function of the Dst index to 
determine the first degree/order spherical harmonic description of the 
near-earth external field and its corresponding induced field. The 
analysis was done separately for data from dawn and dusk with the follow- 
ing results: 

Dusk external : - 20.3 - 0,68 Dst (nT) 

Dusk internal: = 29987.7 + 0.240 q^ (nT) 

Dawn external: q^ - 18.62 - 0.63 Dst (nT) 

Dawn internal: g ]° = 29992.3 + 0.287 qj° (nT) 

Comparison with POGO data indicates that the constant term relating 
q^ 0 and Dst has changed about 20 nT between 1970 and 1980.,. presumably 
due to an increase in the average ring current intensity. 

A local time variation of the external field persists even during 
very quiet magnetic conditions. Both a diurnal and 8-hour period are 
present. A crude estimate of Sq current in the 45° geomagnetic 
latitude. range is obtained for 1966-1970. The current strength, located 
in the ionosphere and induced in the earth, is typical of earlier deter- 
minations from surface data, although its maximum is displaced in local 
time from previous results. 


INTRODUCTION 


Large scale current systems are known to exist on the boundaries of the 
magnetosphere. In the magnetotall and in an equatorial current sheet 
sometimes called the ring current. The fields from these current systems 
are likewise large scale and have been mapped throughout most of the 
magnetosphere and described as functions of the various magnetic activity 
indices and of local time (Fairfield, 1960,1971; Sugiura et al. 1971; Sugiura, 
1972a, b; Sugiura and Poros, 1973, 1979; Hedgecock and Thomas, 1975) The 
equatorial current sheet is apparently the major contributor to these 
magnetic fields at and near the earth's surface. For purposes of this 
paper we cannot distinguish the various sources and will refer to the 
fields of interest as Large Scale External (LSE) fields. These are to 
be distinguished from more localized fields due to concentrations of 
current in "electrojets" or field aligned currents. 

LSE fields are apparently always present. In magnetic observatory 

data their existence is inferred by study of the temporal variations they 

cause in the measured field. When analyzed on a global scale it becomes 

apparent that many of the observatory temporal variations can be described 

in terms of a systematic global morphol ogy (o.g. Sugiura and Chapman, 

1960), particularly during periods of significant magnetic activity. 

Typically, magnetically "quiet" days are selected and used to provide a 

"zero" level or base level from which to measure the temporal variation. 

If we designate the LSE fields by D, then we can write: 

-► -y + 

D = Dst + DS (1) 

-»• 

where Dst is the symmetric portion with respect to the earth's rotational 
axis and DS is the asymmetric port ion. The horizontal component of Dst 


at the equator is the commonly used "Dst index". There is some evidence 
(e.g. Langel and Sweeney, 1971) that OS can be very large during magnetic 
disturbances but tends to zero during quiet times. We doubt, however, 
that it actually goes to zero at any time. These definitions assume that 
DS is magnetospheric in origin, rather than ionospheric (see e.g. Cummings, 1966; 

langel and Sweeney, 1971; Crooker and Siscoe, 1971), and that one can separate 

+ 

DS from the supposedly ionospheric daily variation , Sq. Averages 

of l)S for selected classes of days (e.g. quiet, disturbed etc.) are often 

denoted by SD (see, e.g., Sugiura and Chapman, 1960). 

It is clear that Dst and DS are relative quantities, i.o. their 
values depend upon the quiet level used to define "zero". The absolute 
magnitude of the USE fields at and near the earth's surface has not previously 
been determined. Data from the recent Magsat mission is ideal for such a 
determination. The Magsat spacecraft (Langel et a]., 1982) was launched 
on October 30, 1979, into a twilight, sun-synchronous orbit with inclination 
96.76°, perigee 362 km and apogee 661 km. A Cesium vapor scalar’ and fluxgal.e 
vector magnetometer measured the field magnitude to better than 2 nanotesla (n!) 
and each component to hotter than 6 nl , r.s.s., including orbital and 
attitude determination errors. Magsat remained in orbit until June 11, 

1980, although accurate vector data was obtained only until mid-May, 

1980, Using two days of Magsat data, Langel et al. (1980) derived a 
global spherical harmonic model of the earth's main field, designated 
MGST(6/80), which included a solution for the lowest harmonic component 
of the LSL field. In particular, the external potential determined was 
of the form: 


a(r/a)(q 1 ° P 1 °(0) + [q^cos^ + s 1 1 sin^]P 1 1 (6)} (2) 


where r,0,4> are the standard spherical coordinates, a is the mean radius 
of the earth, and P n m are Schmidt quasi-normal ized Legendre functions. 

The coefficients determined were: s 20.4 nT, q^ 1 = -0.6 nT and 

= -0.4 nT. Previous studies (e.g. Chapman and Price, 1930; Rikitake 
and Sato, 1957; Langel and Sweeney, 1971) have indicated that most of 
the USE field is represented by these first order terms. As noted by 
Langel et al. (1900), the equatorial value of the II component computed Irom 
the first order terms should closely represent the average absolute 
disturbance caused by LSf* currents for the time period of the. data. This 
in turn implies the ability to estimate an absolute baseline for Dst, at that 
time . For the period of data in MGST (6/80 ) the average Dst was 4.6, nT as 
compared with the average absolute value of equatorial H of -20.4 nT. 

The present paper is an extension of the work of Langel et al. 

(1900). Maqsot data has been analyzed for list l)e tween ’ 20 ril and a 
relation between q^ and list established. Because of the Magsat orbit, 
only twilight data is available. This was analyzed separately for dawn 
and dusk to take any local time asymmetry (i.o. 1)5) into account. Extrapolation 
to qj° = 0 then enables us to estimate the true internal field from 
the core and the induced internal field from the time changing external 
fields. The 1 analysis has proceeded in two ways, by spherical harmonic 
modeling of an appropriate Magsat data set and by analysis of individual 
Magsat passes. The two methods yield nearly the same results. 

A re-analysis of magnetic field data from the P0G0 (Polar Orbiting 
Geophysical Observatories, 0G0 -2, -4, and -6, see e.g, Langel, 1974a) 
was carried out for comparison with the results from Magsat. The P0G0 
data yielded less accurate results than the Magsat data since they 


are field magnitude (scalar) data only. Some characteristics of the 
POGO satellites are presented In Table 1. 

Table 1: Summary of POGO Satellite Characteristics. 


SATELLITE 


INCLINATION 


PERIGEE 


APOGEE TIME PERIOD OF DATA USED 


OGO -2 87.30 410 km 1510 km 10/65 - 5/67 
OGO -A 86.00 410 km 910 km 0/67 - 1/69 
OGO -6 82.00 400 kin 1100 km 6/69 - 8/70 


The POGO data also yielded a relation between qj° and Dst. Since this 
data covers all local times, the variation of that relation with local time was 
explored. Comparison with the Magsat results indicates that the zero 
level of Dst has changed about 20 nT between 1970 and 1980. The Sq 
current system is also apparent in the POGO data. 

INITIAL DATA SELECTION FOR SPHERICAL HARMONIC ANALYSIS 

We have begun with the assumption that Dst is a good indicator of 
relative changes in the LSE fields. Dr. M. Sugiura kindly furnished us 
the latest Dst values for the time period of the Magsat mission 
(November, 1979 through June, 1980). These values were derived from 
four observatories (Honolulu, San Juan, Hermanus and Kakioka) and contain 
final baseline adjustments. The Magsat data were then initially screened 
on the three hourly Kp index by choosing only data for Kp < 1- and for 
which the previous three hourly Kp index was <2o. Component data at latitudes 
poleward of 50° geomagnetic latitude was excluded to minimize the effects 



of field aligned and ionospheric currents In the auroral regions, while scalar 
data were retained. After the initial screening, the data was sorted 
into separate subsets for dawn and dusk, according to the Dst 
value at the time of the measurement, and visually scanned both for 
quality and to assure minimization of the effects of short wavelength 
external fields, lable 2 is a tabulation by Dst range of the number of 
data points available for each Dst range and an evaluation of their 
geographical coverage. The geographical coverage from most of these data 
sets is far from satisfactory and, as will be seen, this has an adverse 
effect on some of the analysis. 

SPHERICAL HARMONIC ANALYSIS BY Dst LEVEL 

Twelve of the eighteen subsets of data described in lable ? were 
utilized to derive spherical harmonic expansions of degree and order 
thirteen in constant coefficients, seven in the first time derivative 
coefficients, and one in external terms. Table 3 summarizes- some of 
the results from these models. 

To begin we focused on the relationship between Dst and This 

is plotted, for both dawn and dusk, in Figure 1. The relationship is 
clearly linear over the range of Dst considered. Regression analyses give; 


q,°- 

19.8 - 

* 0.63 Dst, 

nT, 

Dusk 

(3a) 

-O 
1 1 

o 

1! 

18.1 - 

- 0.55 Dst, 

nT, 

Dawn 

(3b) 


with an rms of 0.54 nT and 0.52 nT, respectively. 

The LSE fields vary with time and cause induced currents within the 
earth which, in turn, give rise to internal fields. Presumably the mean 
of such fields, for magnetically quiet times, is absorbed into the g^ 
internal term of standard spherical harmonic analyses of the main field. 


To attempt to Isolate this induced component of internal field from the 
field originating in the core, wo examined the relationship between qj° 
and ijj®, as tabulated in Table 3. figure 2 is a plot of these* variables. 
Clearly the scatter in the determination of qj° is largo. If, however, 
we ignore the three circled points, a trend i$ evident. Regression 
analyses of the remaining points give: 

« -29991.7 + O.Seq^ nT, Dusk (4a) 

gj 0 e -29990.9 + 0.35(1!° nT, Dawn (4b) 

with rms of 1.25 nT and 0.52 nT, respectively. The? circled points are 
for Dst o -10 and 15 nT (Dawn) and Dst » 15 nT (Dusk), from Table 2, two 
of these data distributions are bad and the other is poor, so we feel 
justified in excluding them from the analysis. It is, perhaps, surprising 
that the corresponding q^ 0 values are in such good agreement with the 
other values of Figure 1. One possible explanation is that a poor 
geographic distribution of data implies an even poorer temporal distribution, 
and the secular variation terms are more highly correlated with the 
constant internal terms than with the constant external terms. This 
could result in larger errors in gj° than in q^° for the same model. 

We are not certain how to assign error bars to the data in either 
Figure 1 or Figure 2. The standard errors computed in the least squares 
modeling procedure are, in most cases, less than 0.5 nT„ but we believe 
this underestimates the actual error. As a rule of thumb we feel that, 

fV fi 

for Dst levels with adequate geographic distribution, and g^ 
are probably accurate to '-2 nT. 


COMBINED SPHERICAL HARMONIC ANALYSIS 

In order to verify the results of the previous section and to attempt 
to Improve the accuracy of (3) and (4), the spherical harmonic analysis 
software was modified to permit a direct solution for (3) and (4) using a 
data set selected from ail of the Dst levels. 

The data selection algorithm was applied separately lor the I line 
intervals (1) November - December, 1979, (2) January - February, 19B0, 
and, (3), March - April, 1980, in an attempt to obtain a uniform data 
distribution in both time and space. For each period, and for dawn and 
dusk separately, vector data within • 50° geomagnetic latitude and 
magnitude data poleward of .150° geomagnetic latitude were collected 
into 5° by !>° equiangular bins over the globe. In regions where vectoi 
data were sparse within 150°, available scalar data were retained. 

Within each bin the data was sorted by time and a mean and standard 
deviation calculated. All data with residual (relative to the 0SFC(9/80) 
model, Langel et al. , 1982) greater than 150 nT, and/or greater than 
2u from the mean, is rejected. The desired number of points in each 
bin was selected so as to obtain roughly the same number of points for 
equal area at all latitudes. This was accomplished by specifying a 
maximum of nine values for each vector component in an equatorial b°x5° 
bin, and scaling the number of points in each bin at other latitudes 
by the cosine of the latitude. An equatorial bin value of twenty seven 
was used for scalar data to be retained poleward of jl50°, while six 
was used for scalar data retained within j. 50°. Each of the nine 


Dst intervals was assigned a weight (in nT) ns follows; 

Dst -20 -15 -10 -5 0 5 10 15 20 

Weight 8 8 8 8 8 10 1? 14 18 

If a bin has more data than is desired after applying the above criteria, 

data is rejected first by an Interval skipping algorithm (to maintain 

good temporal distribution) and then by eliminating data with higher 

weights. A computer error resulted in a less than proportionate selection 

from the Dst * -15nT data interval, but the final data sets, dawn and 

dusk, both have adequate temporal and geographic distributions. The mean 

and sigma of the residuals for each data type within each bin for the 

resulting global distribution was computed. To obtain a managable data 

set for least squares modeling while retaining good temporal and geographic 

coverage, the above global data sets were reduced to approximately one 

third size. Within each 5°x5 0 bin. and for each data type, an interval 

skipping algorithm was used (with data sorted by time) to take every 

third point. Any further reduction required was accomplished by eliminating 

that data with the largest absolute deviation from the mean within the 

bin. The data distribution of the resulting global data sets, which are 

used in the combined spherical harmonic analysis, are presented in Table 

4 by Dst interval. 

* 

The spherical harmonic analyses for dawn and dusk using the combined 
data sets from all Dst levels yielded the following results; 


q j° * 18.62 - 0.63 Dst 

nT, 

Dawn 

(6a) 

q x ° = 20.30 - 0.68 Dst 

nT, 

Dusk 

(5b) 

$1° = -29992.3 + 0.29CJ! 0 

nT, 

Dawn 

(Oil ) 

gj° « -29987.7 + O^qj 0 

nT, 

Dusk 

(6b) 


These are to be compared with equations (3) and (4) respectively. We 


consider the agreement to be quite good, and because the data distribution 
Is better for equations (5) and (6) we consider these to be more accurate 
than (3) and (4), The difference In the constant term between (6a) and 
(6b) was unexpected and will be addressed In another paper In preparation, 

ANALYSI S Of INI) I VIDUA! PASS! S 

The residual field after subtraction of a core field model (spherical 

harmonic analysis) from the measured data, o.g* 

> ► ► 

AB - B(moasured) - B(computed) (’/a) 

> * 

aB o ]B (measured)} - | B (computed ) | (7b) 

Is due to (i) fields from external currents and the corresponding induced 
currents, (ii) crustal anomalies and, (in), measurement and field model 
errors. We will assume that (ill) is negligible. At low latitudes, away 
from noon, the external portion of (7) is dominated by the long-wavelength 
L5f; fields. Wo also assume that (ii) is negligible at these wavelengths, 
and that (1) is described by a potential field of the form (Chapman and 
Price, 1930; Rikitake and Sato, 1957): 

V » a [(r/a)e + (a/r ) 2 1 )]cosu (8) 

where r is •geocentric radius, o is geomagnetic, latitude and "a" is the 
mean radius of the earth. Here e corresponds to in {?,)> except 
one uses geographic and the other geomagnetic latitude. Lange 1 and 
Sweeney (1971) and Langel (1974b) showed that a pass of satellite data 
between ±45° latitude was sufficient to solve for e and i in a 
least squares sense, even if only scalar data is available. Such a 
solution is more accurate when component data is available, as from 
Magsat. 


We have computed e and i for all Magsat passes, both from the vector 
and from the scalar data. The residual fields used were relative to the 
MOST ( 4/81 ) spherical harmonic model , This model is of degree and order 
thirteen in constant terms and seven In first time derivative terms, 
Though unpublished, MGST(4/81) may be obtained from the authors. 


Figures 3 and 4 show e plotted versus Dst for dawn and dusk, The 
range of Dst is limited to ±20 nT; there are 1211 points at dawn and 1 16 i 
points (. dust, Uoijross t on lines fit In I he dale give: 


e z 21.3 - 0,62 Dst 

nT, Dawn 

(0a) 

(' 22.4 - 1). /? Dst 

fli, Dusk 

('«’) 

with standard deviations of 3.9 

nT and 4.2 riT, 

rospocti voly, and corrolnf I 

coefficients of 0.82 and 0,85, 

respecti vely. 


Figures 5 and 6 show 1 plotted verses e. 

Regression linos fit to 

the data give: 



J « -6,1 + 0. I7e 

nT , Dawn 

(10a) 

i a 0.3 4 0.21e 

nT, Dusk 

(10b) 


with standard deviations of 3.6 nT and 3.8 nT and correlation coefficients 
of 0.31 and 0.39, respectively. 

Equations (9) are in very good agreement: with equations (3) and (5) 
and so we 4re confident we have established a clear relationship between 
q}° and Dst for the epoch of the Magsat data. 

The relationship between i and e is less satisfactory, although the 
slopes in (10) are certainly consistent with those in (4) and (6). We had 
expected the correlation coefficients to be larger and the scatter iri 
Figures 5 and 6 to be less. Theoretically, the i/e ratio will depend 
upon the rate of change of e and we thought that taking this into account 
might improve the results. However, a plot of i/e as a function of the 


rate of change of e showed no correlation. The computed correlation 
coefficient was 0.06. 

Another difficulty with (10) arises from the field model, MG5T(4/81), 
used to reduce the data. No account was taken of the external or Induced 
fields in deriving this model. The effect of this on the determination 
of e and 1 is not known, but clearly the term in the MOST (4/81) 

model will have absorbed some of the induced fields so that the constant 
terms in (ID) are somewhat uncertain. 


ANALYSIS OF POGO DATA 

The POGO data are scalar field measurements only. Attempts to 
solve for external fields with spherical harmonic analysis using this 
data have been inconclusive. However* analysis of Individual passes to 
solve for the coefficients e and i in (10) is possible (Langel, 1974b). 

To carry out this analysis the POGO data were subdivided into Dst bins 
5nT wide centered at 5 nT intervals between '20 nT, into local time 
bins 2 hours wide centered at 0, 2, 4, ...22 hours, and into geomagnetic 
latitude bins 5° wide centered at even 5° intervals between 1 40°. 

The residuals relative to the P0G0(8/71) spherical harmonic model {Langel, 
1974a) were averaged within each bin. Solutions for e and i were then 
obtained for each Dst and local time. For each local time a regression 
analysis of the form: 

e - a -*• b Dst (11) 

and 

i » c + d e (12) 

was computed. Table 5 gives the values of a, b, c and d for each local 
time; Figures 7 and 8 contain the corresponding plots. Notice that the 
zero hour value is present on both ends of the plot. A local time 
variation is present in a, b and c, but we regard d to be constant at 
its average value, 0.22, to within the accuracy of our determination. 

The local time dependance of a and b is clearly of a sinusoidal 
nature. Fourier analysis yields: 

a = 0.40 - 2.87sin(t~8.2°) + 1.82sin(3t+12.05O), nT (13) 
b - -0.600 + 0.134sin(t*il3°) (14) 

where t is the local time in degrees. Calculated values from equations 
13 and 14 are given as solid lines in Figure 7. 


A diurnal variation of DS was found by Sugiura and Chapman (I960), 
who also found significant semi-diurnal variation. To our knowledge the 
period of 8 hours has not been seen previously in fields originating 
external to the ionosphere. Although small, we see no reason at present 
to doubt its reality. The absence of the 8 hour period in b means that 
its magnitude is essentially constant regardless of Dst level, for 
Dst within 20 nT. This seems to suggest that its origin lies in 
magnetopause currents rather than the ring current. 

The semi-diurnal terms in the Fourier analysis were; 
a (semi-diurnal) - 0.72sin(2t + 53°) (15a) 

b (semi-diurnal) = 0.23sin(2t + 105°) (15b) 

However the spectral power density of the semi-diurnal term for a is a 
factor of seven down from the 8 hour period and for b it is a factor of 
35 down from the diurnal term. We do not regard these terms as significant, 
although their amplitudes, relative to the diurnal coefficients in (13) and 
(14), are comparable to those found by Sugiura and Chapman (1960). 

It is apparent that the local time variations persist at all Dst 
levels and that the strength of the diurnal variation, in particular, 
grows as Dst becomes more negative. It. should be kept in mind that the 
present a/iaylsis considers data averaged over Dst. The results, therefore, 
indicate the variation fields averaged over both Dst and time and are 
thus unable to distinguish the variations of a and b over short periods 
of time. That such variations exist has been shown clearly by Crooker 
(1972). 


The potential 


V = a (a/r)2 i cose (16) 

describes fields originating internal to the satellite altitude. Such 
fields may be of two origins: from currents in the ionosphere or from 
currents induced in the earth. If 

i * i in + i 5 » (17) 

where i^ represents fields induced by currents in the magnetosphere and 
i s represents ionospheric currents and their corresponding induced earth 
currents, (presumably Sq for these Dst levels), then: 
im + is “ c + d ©• (18) 

If the external field, e, goes to zero 1 m should also go to zero from 
which we infer that c is strictly associated with Sq* It would be nice 
if we could say that d is strictly associated with earth currents induced 
by inagnetospheric currents, as is perhaps likely, but that cannot he 
concluded with certainty. Any portion of Sq which varies with Dst will 
also be included in d. 

Under the assumption that contributions of Sq to d are small, we can 
substitute c for i in equation (16) to obtain the potential for Sq. 

There is a zero-level uncertainty to this procedure because some portion 
of the Sq field will have been absorbed into the spherical harmonic 
model used to represent the core field. Table 5 gives the coefficients 
of this potential for each of twelve local times. These potentials are 
approximations since they are derived using data only between i: 4 0 0 geo- 
magnetic latitude along a single local time meridian, and assume a strictly 
P 1 °(cose) functional relationship. This procedure does not yield a 


global potential. It is possible, however, to compute a current function, 
along each of the local time meridians, between 45° latitude. 

Under the assumption of zero current near midnight, we have added 2.0 to 
the values of Table 5 when performing this calculation. 

The resulting current is shown in Figure 9. Coordinates are 
geomagnetic latitude and solar local time. Comparison with typical Sq 
current systems derived from observatory data indicates good agreement in 
current direction and magnitude. Those systems, however, indicate that 
the current maximum should occur near 1100 hours wheras the maximum 
current in Figure 9 occurs near 1300 hours. We are unable to account for 
the discrepancy at this time. 

Mead (1964) and 01son(e.g. 1970) have pointed out that the LSF. fields 
at the earth contain variations similar to the Sq field. This has raised 
the question , still unresolved, as to how much of the Sq field originates 
in ionospheric currents and how much originates external to the ionosphere. 
The equivalent current system of Figure 9 must reside in the ionosphere 
and the earth itself, at least to within the accuracy of our determination. 
This would indicate that the major part of Sq is ionospheric in origin, 
although it does not exclude a magnetospheric contribution. 


COMPARISON OF MAGSAT AND POGO RESULTS 
Writing 

~ a + b Dst (19a) 

i » d q^, (19b) 

a summary of results for dawn and dusk is given in Table 6. Comparing 
Magsat and POGO shows that the asymmetry between dawn and dusk is present 
in both data sets. The values of d from POGO are identical with the 
Magsat values from the combined spherical harmonic analysis, which we 
regard as the most accurate Magsat determination. Variation of b is 

larger in the POGO results than in Magsat and the values of a are 

significantly different. 

The most important difference is that in the value of a. Magsat 

data yields a value 23.1 n t higher than POGO data at dawn and lb, 6 ni 

higher at dusk. Thus, from (19a), for the same level of Dst the 

magnitude of the external field near the earth is about 20 nT higher at 

1980 than during 1966-1970. This difference seems somewhat excessive. 

The question arose as to whether use of the P0G0(8/71) spherical harmonic 

model to reduce the POGO data might be introducing a bias since it is 

only of degree/order ten. Accordingly, a portion of the POGO analysis 
* 

was redone using the GSFC(9/80) model (Langel et al., 1982). No sig- 
nificant change in the a or b coefficients was found. The data in both 
these models for the 1966-1970 time period is predominantly scalar 
(from POGO), which still leaves some question as to how accurate a model 
is obtainable. However, we have no indication that the models could 
be in error by as much as 20 nT, 

Assuming the calculation to be correct, we are led to conclude that 


either the definition of the Dst baseline has changed or that the sources 
of the external field have changed. The baseline for Dst is derived 
from selected magnetically quiet days, after removing effects from Sq. 

The method of analysis has been self consistent for a period of time 
exceeding that under discussion here (Sugiura, private communication). 

We see no reason to attribute any change to its definition. On the 
other hand, magnetic quiet is defined in terms of the temporal variation 
of the field at the various observatories, and a long-term change in 
external fields might not be detectable as such in surface data but 
could very well appear to be part of the main field secular variation. 

Let us assume that the external sources have changed. As seen from 
near the earth these external fields may be considered to be the result 
of two competing mechanisms. The ring current and tail current, or 
equatorial current sheet, inflates the magnetosphere. An increase in 
the average level of these currents would result 'in an increase in the 
value of the coefficient a. Compression of the magnetosphere by increased 
solar wind pressure (an increase in magnetopause currents) would result- 
in a decrease of the coefficient a. Between 1966-1970 (POGO epoch) and 
1980 (Magsat epoch) we observe the coefficient a increasing by about 20 
nT which tmplies increased inflation of the magnetosphere or decreased 
magnetopause compression. In fact, on average, the magnetopause compression 
increased. From Mead (1964) we find that the field at the earth from 
Magnetopause compression increases by about 7.5 nT for each earth radius 
of contraction of the radius of the subsolar magnetopause. It is commonly 
accepted that the subsolar magnetopause varies inversly with the sixth 


power of solar wind pressure. The average solar wind pressure for the 
POGO epoch was 1.54 x 10“ 8 dynes/cc and for the Magsat epoch 2.3 x 
10~ 8 dynes/cc (J. King, private communication). Thus the magnetopause 
radius, on average, was less during the Magsat period than the POGO 
period. A simple calculation shows that the field at the earth would have 
decreased (or the coefficient, a, would have increased) by a little over 
5 nT. 

This implies that the amount of inflation of the magnetosphere due 
to the ring and tail currents was greater in 1900 than in 1966-1970. 

The four year period of the POGO data is considerably longer than 
the 7.5 months of Magsat data. To determine if the coefficients a and b 
changed during 1966-1970, the POGO data was subdivided into nine time 
intervals and separate analysis made. No apparent change in the coefficients 
was found. 


SUMMARY AND DISCUSSION 


Estimates of the magnitude and temporal variation of the near-earth 
I.SF fields have been available for some time. Magnetospherie models have 
been designed to represent the magnetic field throughout the entire 
magnetosphere, of which the vicinity of the earth is a very small part. 
Further, the fields of interest in those models are very weak near the 
earth so that even small model shortcomings can result in a large percentage 
error in representing the near-earth field (Olson and Pfitzer, 1981), 
Complicating the issue further is the difficulty inherent in attempting 
to isolate the external fields from the internal fields near the earth. 
Observatory data has been used successfully to show that most of the Sg 
fields are external to the earth and the most plausible theories regarding 
origin of the Sq field indicated currents in the ionosphere. Now, however, 
we find (Olson, 1970) that a portion of Sq may ho niagnetosphoric in 
origin. Plausible arguments were made from observatory data that, at 
least during magnetic storms, the symmetric part of I) (Dst) originated 
in a ring current external to the ionosphere. IK fields were originally 
analyzed iri terms of ionospheric currents (e.g. Chapman, 1939). later 
it was argued that at least some of DS was magnetospherie rather than 
ionospheric (e.g. Akasofu arid Chapman, 196d). However definitive 
separation between ionospheric and magnetospherie sources requires data 
at altitudes between the two, i.e. from a satellite (e.g. Lange.] and 
Sweeney, 1971) and, where possible, data in the actual current, region 
(e.g. Sugiura, 1972b). Study of LSE fields with observatory data has 
focused on the time-varying portion of those fields, especially during 
magnetic storms. From measurements at the earth's surface it has not 


been possible to describe LSE fields which may be nearly constant In 
time. Field variations which correlate with solar activity have been 
found(Yukutake, 1965; Alldredge, 1976; Suglura, 1976; Harwood and Malln, 
1977; Suglura and Poros, 1977; A1 1 dredge et a 1 . , 1979; Yukutake and 
Cain, 1979) and attributed to external fields, although their Isolation 
Is extremely difficult and the accuracy with which the determination is 
made is problematical. 

We have derived a description of the average LSE field and the 
corresponding Induced field between 45° geomagnetic latitude for 
magnetically quiet periods in terms of a first degree/order spherical 
harmonic analysis of Magsat data. This model relates the LSI; field 
directly with the Ost index, establishing an absolute baseline for Dst at 
1980. Separate analyses for dawn and dusk local times showed a distinct 
asymmetry. This was confirmed by an analysis of POGO data from 1966-1970 
in which the asymmetry was found to have both an 8 and a 24 hour variation 
with local time. The amplitude of the local time variation was greater 
in the POGO data than in the Magsat data, even though the Ost levels 
were the same. Both sets of data show that the amplitude of the asymmetry 
Increases as list becomes more negative. Thus, although the external 
field asymmetry is greater during magnetic storms, it never entirely 
disappears. We are unable to tell if this is because the ring-current 
never becomes completely symmetric or if the source of the asymmetry 
lies in magnetopause and magnetotail currents. 

Previously Langel et al, (1980) found q^ 0 to be 20.4 nT for provisional 
Dst = 4.6 nT. The adjusted Dst for the epoch of that analysis is 3 nT. 

From (5) we now predict q^ to be 16,7 nT at dawn and 18.3 nT at dusk, or 
about 2-4 nT lower than the earlier result. The statement by Langel et 


al. (1980), "To estimate the longitudinally symmetric portion of the 
disturbance one should subtract 25 nT from the value of Dst..." is 
Incorrect* Equations (5) and (6) give the correct estimate. 

Comparison of POGO and Magsat results indicates an Increase In the 
average quiet-time intensity of the LSE fields from 1966-1970 to 1980, 
which must be attributed to an intensification of the ring current. The 
magnitude of the change is about 20 nT. Because the POGO analysis is 
dependant on scalar data only, it is difficult to estimate the accuracy 
of this result. However, comparison of similar calculations from Magsat 
showed that separation of the internal and external fields on individual 
passes using only scalar data was usually within 4 nT of the separation 
using only vector data. 

Sq fields have also been shown to he discernable in the POGO data. 

To the accuracy of our determination, these arise entirely internal to 
the spacecraft altitude, suggesting that any contributions to Sq from LSf 
fields is small. 

finally, the magnitude of the coefficient of Pj° for induced fields 
is about 0.24 to 0.29 of the corresponding coefficient for the inducing 
field. This has implications for the investigation of upper mantle 
conductivity as already studied for the POGO data by Langel (1974b) and, 
more extensively, by Pidwall (1981). 


TABLE 2; DATA DISTRIBUTION FOR EACH DST INTERVAL 
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ABLE 3: PARAMETERS FROM SPHERICAL HARMONIC MODELS BY Dst, ALL UNITS ARE nT 


*ir 


TJ 

0) 

U 

11 

X3 

e 

*r» 

(/) 

1 

r ■' 

Cl 

<0 
I » 

o 

u 


ra 


c 

Ot 

00 

oo 

at 

io 

f-* 

at 

CM 

O 

rH 

oo 

Cw 

• 

• 

• 

• 

« 

• 

♦ 

♦ 

» 

t 

* 

<u 

00 

m 

io 

«d* 

CM 

t-l 

r-l 

oo 

CM 

CM 

at 

+j 

U3 

60 

CO 

co 

CO 

r*» 


60 

00 

00 

r**. 

C 

at 

at 

at 

at 

at 

at 

at 

at 

at 

at 

at 

*r* 

at 

at 

at 

at 

at 

at 

at 

at 

at 

at 

at 


CM 

CM 

CM 

CM 

CM 

CM 

; i 

CM 

CM 

CM 

CM 


rfl 
tjil ' 


2= 

< 

a 


fO 

c 

L. 

Of 

+J 

x 

OJ 


O' 


«V 

ST) 

O 


H 

CM 

rH 

cn 

at 

rH 

m 

CO 

M‘ 

CM 

CO 

• 

• 

• 

* 

• 

• 

• 

* 

• 

V 

• 

o 

CM 

LO 

hn 

rH 

ro 

at 

Kt 


O 

rn 

pH 

rH 

rH 

rH 

<N 

CM 


y*H 

pH 

CM 

CM 


+J 

</> 

Q 


IT) 


o 


LO 


o 


in 

( 


in 

rH 


O 

t H 


m 


o 


in 

i 


-10 25.5 -29975. 


TABLE 4: DISTRIBUTION OF DATA BY Dst INTERVAL IN COMBINED DATA SETS 
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TABLE 6: SUMMARY OF RESULTS 
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Captions 

1: First degree/order coefficient of the external potential, q^, 
as a function of the Dst index for dawn and dusk local times, 
Each coefficient qj was determined in a spherical harmonic 
analysis of the main and external fields in which all data 
was from periods of time when Dst was within 2.5 nT of its 
nominal value. 

2: First degree/order coefficient of the internal coefficient, 

Oj , as a function of the first degree/order external coef- 
ficient, . Each determination uses data for which Dst 
was confined to a *2.5 nT range, as described in the caption 
to figure 1. 

3: First degree/order coefficient, e, of the external potential 
fit to the X and Z component residual fields from individual 
Magsat passes verses the Dst index at the time of the equator 
crossing of that pass. Data are for dawn local time only. 

4: Same as Figure 3 except dusk local time only. 

5: First degree/order coefficient, i, of the internal field 
versus first degree/order coefficient, e, of the external 
field as determined from the X and Z component residual 
fields from individual Magsat passes. Data are for dawn 
loc.il time only. 

6: Same as Figure 5 except dusk local time only. 

7: Local time variation of coefficients a and b from equation (11) 
describing the external field for POGO data as a function 
of Dst. 

8: Local time variation of coefficients c and d from equation (1?) 
describing the internal disturbance field potential for POGO 
data as a function of Dst. 

9: ‘Equivalent current system corresponding to the disturbance 
field internal to the POGO orbit and extrapolated to the 
condition of zero field from magnetospheric sources. 

2.5x10^ Amps, flows between streamlines. 
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